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px| ABSTRACT 

I I Context. The CoRoT space mission has been searching for transiting planets since the end of December 2006. It has led to the detection 

of about twenty Jupiter-size planets and three planets with radius < 5 R^. The latter are CoRoT-7b, the first super-Earth observed 
in transit, and two validated Neptunes, CoRoT-24b and c, in a multiple system. 

Aims. We aim to investigate the capability of CoRoT to detect small-size transiting planets in short-period orbits, and to compare the 
number of CoRoT planets with 2.0 < Rp < 4.0 Re with the occurrence rate of small-size planets provided by the distribution of Kepler 
planetary candidates (Howard et al. 2012, ApJS, 201, 15). 
^ Methods. We performed a test that simulates transits of super-Earths and Neptunes in real CoRoT light curves of six long observational 

^ runs and searches for them blindly by using the transit detection pipeline developed at the Laboratoire d'Astrophysique de Marseille. 

Results. The CoRoT detection rate of planets with radius between 2 and 4 Rg, and orbital period P < 20 days is 59% (31%) around 
stars brighter than r' = 14.0 (15.5). The vast majority of the missed planets went undetected because of a low transit signal-to-noise 
I I ratio (5 /N). However, in some cases, additional instrumental or astrophysical noise may prevent even transits with relatively high 

S /N, i.e. S /N > 10, from being revealed. By properly taking the CoRoT detection rate for Neptune-size planets (2 < Rp < 4 R^) 
and the transit probability into account, we found that according to the Kepler planet occurrence rate, CoRoT should have discovered 
^ 12 ± 2 Neptunes orbiting G and K dwarfs with P < 17 days in six observational runs. This estimate must be compared with the 

1/^ validated Neptune CoRoT-24b and five CoRoT planetary candidates in the considered range of planetary radii, the nature of which 

is still unsolved. We thus found a disagreement with expectations from Kepler at 3 cr or 5 cr, assuming a blend fraction of 0% (six 
00 Neptunes) and 100% (one Neptune) for these candidates. 

Conclusions. This underabundance of CoRoT Neptunes with respect to Kepler may be due to a) an underestimate of the uncertainty 
on the Kepler planet occurrence; b) an underestimate of the false-positive probability of the Kepler small-size planetary candidates; 
c) an overestimate of our prediction of CoRoT Neptunes related to the number of G and K dwarfs observed by CoRoT; or d) different 
stellar populations probed by the two space missions. Regardless of the origin of the disagreement, which needs to be investigated in 
more detail, the noticeable deficiency of CoRoT Neptunes at short orbital periods seems to indirectly support the general trend found 
in Kepler data, i.e. that the frequency of small-size planets increases with increasing orbital periods and decreasing planet radii. 

Key words, planetary systems - techniques: photometric - stars: solar-type - stars: statistics 

1. Introduction About a hundred low-mass planets have been recently dis- 

covered by the HARPS and HIRES radial-velocity surveys (e.g.. 



More than 750 exoplanets have been found since the discov- [Mayor et al. [201 Ijpbward et al. |2010[ and references therein) 



ery of the extrasolar planet 51 Peg b around a solar-like star by From these early discoveries, low-mass planets seem to be abun- 
[Mayor & Queloz ( 1995 ). Most of them are gaseous giant planets, dant. Focusing on 376 solar-type stars (from late-F to late-K 



but much effort is currendy put into detecting small-size planets, dwarfs ) observed with the HARPS spectrograph, [Mayor et aT 
such as Neptunes, super-Earths, and Earth-sized planets, with (j201 1) found that about 28±5% of them host a planet with a min- 
radial-velocity and space-based transit surveys. imum mass Msin / = 3 to 30 M® and orbital period f < 50 days. 

With a more limited sample of G and K dwarfs observed with 
the HIRES spectrograph, [Howard et al. |p010| ) reported an oc- 
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currence rate of 18 + 5%, slightly lower than the rate of |Mayor] 
et al. | (201 1 1, for the same mass domain and period range. 

The observed occurrence rate of low-mass planets around 
dwarf stars, its dependence on the properties of the host stars 
(metallicity, stellar mass, etc.) and on the planetary parame- 
ters (mass, semi-major axis, and eccentricity) represent a pow- 
erful test for models of planet formation and migration (e.g., 
Mordasini et al. ||2009al IMordasini et al. „2009b Mordasini et| 



al. [2012 and references therein; Mayor et al. |201 1 1. 

In addition to their contribution to general statistics, transit- 
ing low-mass planets are extremely interesting to study because 
one can measure their radius and mass free of the sin / degen- 
eracy if radial-velocity variations induced by the planet on its 
parent star or transit timing variations in a multiple system are 
observed. This enables a detailed characterisation of small-size 
planets and studies of their internal structure. 

The transits of small-size planets can be observed from the 
ground only if their host stars have sizes considerably smaller 
than that of the Sun. This is, for example, the case of the super- 



Earth or mini-Neptune GJ 1214b (Charbonneau et al. 2009 1. 



In general, we have to move to space to detect small-size plan- 
ets around solar-like stars as the photometric precision from the 
ground is insufficient. 

1.1. The CoRoT space mission. 

The CoRoT space mission is the pioneer for the detection of 
transiting planets from space. It is operated by CNES[^ with 
contributions from Austria, Belgium, Brazil, ESA, Germany 
and Spain (Baglin 2003). CoRoT was launched at the end of 
December 2006 on a polar orbit at an altitude of ~ 900 km. 
With its 27 cm telescope, it has been monitoring the optical flux 
of thousands of stars with 11 < V < 17 for ~ 150 days at the 
longest, in several fields of view pointing mainly in two different 
directions: one towards the Aquila and the other one towards the 
Monoceros constellation. The field of view of the telescope for 
the exoplanet science was ~ 4 deg^ before one DPU on board 
was lost through technical problems on 8 March 2009. Since 
then, only one CCD dedicated to exoplanet science is currently 
working, giving a halved field of view of 1.38 deg^. We refer the 
reader to Auvergne et al. (2009 1 for a detailed description of the 
instrument. 

To date, CoRoT has led to the detection of three small-size 
planets with Rp < 5 R(s : CoRoT-7b, the first super-Earth ob- 



served in transit (Leger et al. 2009 1, and two transiting Neptunes, 
CoRoT-24b and CoRoT-24c, in a multiplanet system ( Alonso et 



I al. [su bmitted). For the last two planets, radial-velocity observa- 
tions performed with the HARPS spectrograph did not have suf- 
ficient precision to derive their masses, given that their parent 
star is faint. However, by using a procedure similar to that de- 
scribed in |Torres et al.| ( |201 l|l, these candidates hav e been "vali- 
dated" as low-mass planets ( Alonso et al. |submitted| l. This means 
that the probability that the observed transits are due to astro- 
physical false positives, such as blended eclipsing binaries, is 
very low. In other words, the planet scenario is most likely true 
( [Torres et al.T20TT] i. 



156, 000 stars with 9 < V < 16 in a fi xed stellar field of 1 15 deg-^ 
located in the Cygnus constellation (j Borucki et al.|[2006 Koch 
et al. [20 To I. With its 95 cm telescope and 42 back-illuminated 
CCDs, it reaches a photometric precision of ~ 50 ppm on a 
my = 1 1 .5 star for a ~ 30 min integration time, while the corre- 
sponding RMS of CoRoT is ~ 200 ppm. Moreover, since Kepler 
has been observing always the same stellar field, as time goes by, 
it has the capability of detecting shallower and shallower transits 
because the signal-to-noise ratio (S /N) scales as the square root 
of the number of transits. Thanks to the unprecedented photo- 
metric precision of the instrument, the Kepler team found 2321 



planetary candidates in the first sixteen months of data (Batalha 
et al. 2012) . The vast majority of them have radii below the 
Neptune radius. 

Up to now, the Kepler team has announced the discovery of 
42 planets with radius smaller than 5 7?®, many of which belong 
to multiple system^ We recall the exceptional discoveries of 
the Kepler- 11 and Kepler-20 systems with six and five small- 
size planets, respectively. In many cases, for instance Kepler-20, 
the planets have been only "validated" because their mass has 
not been measured either through radial-velocity observations 
or transit timing variations (e.g., Tor res et al.||2011 [Fressin eT 
al. [201 la ). Kepler-20e wi th Rj, - 0.87 is the smallest planet 
found by the Kepler team ( Fressin et al. |20lTb| . 

From stellar population synthesis and Galactic structure 
models, Morton & Johnson] ( |201 1) claimed that more than 90% 
of Kepler planetary candidates are very likely real planets. Based 
on this fiducial lower limit, Howard et al. (2012) have per- 
formed detailed statistical studies on the Kepler planetary can- 
didates. After correcting for the alignment probability, these au- 
thors derived an occurrence rate of 13 ± 1% for planets with 
radius 2 < Rp < 4 7?®, orbiting G and K dwarfs in P < 50 days. 
Moreover, Howard and coworkers found that the planet occur- 
rence increases with decreasing planet radius and increasing or- 
bital period. 

1.3. Outline and purpose of tlie present work 

In the present work, we first of all aim to investigate the ca- 
pability of CoRoT of detecting small-size planetary transits. To 
that purpose, we simulated planetary transits of super-Earths and 
Neptunes in real CoRoT light curves of G and K dwarfs and 
searched for them blindly, that is, without knowing in advance 
in which light curves the transits had been simulated (Sect. |2]i. 
We used the transit detection pipeline currently operating at the 
Laboratoire d'Astrophysique de Marseille (Sect.[3]l, which is in 
charge of the CoRoT alarm mode. This allowed us to derive the 
CoRoT detection rate for planets with different sizes, from 1.3 
to 5.0 R@ (Sect. 4.1 1, and to investigate the obstacles to the de- 
tection of small-size planets with CoRoT (Sect. |4.2| i. By deriv- 
ing the CoRoT detection rate of planets with radii from 2.0 to 
4.0 R^, we compared the number of small-size planets discov- 
ered by CoRoT with expectations from the Kepler occurrence 
rate provided by 'Howard et al. ( 2012) l, after taking the transit 
probability into account (Sect. 4.3 1. 



1.2. Tlie Kepler space mission. 

The other space-based mission aimed at detecting planetary tran- 
sits is the NASA Kepler mission, launched in April 2009 in 
a heliocentric orbit. It performs photometric observations of ~ 

' French space agency: Centre National d'Etudes Spatiales 



2. Simulations 

We made use of real light curves (N2 level data) of the CoRoT 
long runs lasting more than 110 days, i.e. LRaOl ( [Carone et al. 
. 2012 ), LRa02, LRa03 ( |Cavarroc et J~||20T2| , LRa04, LRcOl 

^ http://kepler.nasa.gov/Mission/discoveries/ 
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( [Cabrera et ai~||2009| l, and LRc02 f\ They refer to six differ- 
ent fields of view observed by CoRoT. For all runs except for 
LRa03 and LRa04, both of the two CCDs dedicated to exoplanet 
science were operating. Out of the ~ 55 000 stars observed in 
the above-mentioned long runs, almost 14000 were classified as 
dwarfs with spectral type G and K, according to the photomet- 
ric classification of Exodat ( jPeleuil et al.|2009| . The remaining 
targets were ranked as either evolved or hotter stars. Among the 
~ 14000 G and K dwarfs, only those brighter than r' - 15.5 
were selected, that is, 9526 stars. From this sample, light curves 
of binaries or planetary candidates found by the CoRoT detec- 
tion team were excluded. This yields 9333 light curves for our 
blind test. 

Here we must point out that from a comparison between the 
Exodat photometric classification and a massive spectral anal- 
ysis of more than 1 000 CoRoT stars wi t h FLA MES/GIRAFFE 
multi-fiber observations, [Gazzano et al.| ( |2010[ ) found that the 
former tends to systematically underestimate the number of 
dwarfs, while there is a fairly good agreement between the spec- 
tral types determined by the two methods. Figure 8 in Gazzano et 
al. ( 2010 1 clearly shows this behaviour for the LRaOl and LRcOl 
stellar fields: the spectral analysis found 52% of the dwarfs while 
only 32% of the stars were classified as such by the photometric 
classification, which is a difference of a factor 1.7. The num- 
ber of subgiants identified by both methods is similar, while the 
amount of giants was clearly overestimated by Exodat. However, 
since the classification based on spectroscopy is only available 
for about a thousand targets in the LRaOl and LRcOl, we could 
only rely on the Exodat classification for the light curve selec- 
tion. Nevertheless, we will consider a correction factor to take 
the underestimate of dwarf stars by the Exodat classification into 
account in Sect. 14.31 

The choice of restricting ourselves to the light curves of G 
and K dwarfs allows us to compare the rate of small-size plan- 
ets detected by CoRoT with the statistics provided by the dis- 
tribution of Kepler planetary candidates. In addition, G and K 
main-sequence stars with r' < 15.5 are the most interesting tar- 
gets for the detection of small-size planets. The fainter the stars 
and/or the bigger their radius, the more diflicult the discovery of 
Neptune s and super-Earths. 

Simulated transits of small-size planets, super-Earths and 
Neptunes, with radius Rp ranging between 1.3 and 5.0 Rg^ and 
orbital period P from 0.7 up to 20 days were inserted into the se- 
lected 9333 CoRoT light curves. A flat distribution was adopted 
both for the planetary radius and the orbital period. To obtain a 
sufficiently large statistics for the detection of super-Earths, we 
repeated our experiment twice. First we simulated only transits 
of planets with 1.3 < Rp < 2.5 7?® in the 9333 light curves. Then, 
we inserted transits of bigger planets with 2.5 < Rp < 5.0/?® into 
the same original light curves. 

We used the universal transit modeler (Deeg 2009 ) to simu- 
late the planetary transits assuming, for simplicity, circular plan- 
etary orbits, the stellar radius and mass of a G2V star for the G 
dwarfs, i.e. Ri, = 1 Rq and M^, - 1 Mq, and those of a K2V 
star for the K ones: R^ = 0.8 Re and = 0.7 Mq (Cox^gOOO}. 
A linear limb-darkening law was adopted with the linear coeffi- 
cient for the CoRoT bandpass u - 0.62 or 0.69 for the G or K 
spectral types, respectively (Sing 2010). The impact parameter 
was chosen at random between and 1 to take the distribution 
of the orbital inclination into account. The initial epoch of the 
transits was also chosen randomly. 



Our simulations do not take planets around one component 
of a binary or triple physical system into account. Such planets 
would be more difficult to detect because their transit depth is di- 
luted by the stellar companion(s). However, for the time being, it 
is not clear how the frequency of small-size planets is affected by 
the multiplicity of stellar systems and as long as we lack this in- 
formation, we will consider only single stars. Moreover, even the 
distribution of the planetary radii of Kepler candidates, which 
was used by Howard et al. | ( 2012| l to determine the planet oc- 
currence as a function of orbital period, implicitly assumes that 
they orbit single stars ( [Borucki et al. |20fT [Batalha et al. 2012) l. 

Another possible cause of dilution of the transit depth is the 
contamination of background stars inside the CoRoT aperture 
mask of a star ( Deeg et al.|2009| l. However, this can be readily 
neglected for our study, because this contamination is typically 
only ~ 5% of the flux collected by CoRoT for a given star (e.g., 
IBonomo et alpOIOl ). 



3. Data analysis 

3. 1 . The LAM transit detection pipeline 

Transits were searched for by means of the transit detec- 
tion pipeline developed at the Laboratoire d'Astrophysique de 
Marseille (LAM). The CoRoT detection team uses seven differ- 
ent techniques to detrend the light curves and search for transits. 
Blind tests organized within the CoRoT consortium to compare 
their performance on real CoRoT light curves have shown that 
the LAM pipeline and two other algorithms have the best perfor- 
mance and a similar detection rate (R Barge, priv. comm.). Our 
pipeline consists of eight steps: 

a) Outliers caused by proton impacts during the passage of 
the satellite across the South Atlantic Anomaly ( jAuvergne et al.| 
2009) are filtered out with a 5-sigma clipping. 

b) Data points oversampled at 32 s are re-binned at 512 s. 

c) Low-frequency variations (i.e., stellar variability) are re- 
moved with a high-pass filter, consisting of a sliding median fil- 
ter with a window extension of 0.5 days (see Sect. |3.3| for a dis- 
cussion on the choice of the filtering time scale). 

d) High-frequency variations are removed with a low-pass 
filter, consisting of a Savitzky-Golay filter ( [Press et al.||1992| 
with a time scale of ~ 1 hour Most of these high-frequency 
variations are related to the signal at the satellite orbit of 6184 s 
(1.72 h). 

e) Jumps/discontinuities that are mainly caused by hot pix- 
els ( Srour et aL][2003 cf . 'Auver gne et al.[[2009 for the specific 
case of CoRoT) or pointing displacements are automatically de- 
tected. After that, data points that fall 0.4 days before and after 
the time of a given hot pixel are removed. These discontinuities 
in CoRoT light curves may represent a major concern for de- 
tecting shallow planetary transits (see below Sect. 4.2 1. Indeed, 



data available at http://idoc-corot.ias.u-psud.fr/ 



after the filtering process (steps a-b-c-d), they give rise to dips in 
the residuals that can be erroneously identified as transits by the 
detection algorithms, such as the box fitting least squares (BLS, 
[Kovacs et al.T2002| . 

The detection of hot pixels is carried out by using a slid- 
ing window of eight data points sampled at 512 s and comput- 
ing the standard error of the mean inside the moving window. 
When a discontinuity takes place in the light curve, the stan- 
dard error suddenly increases, allowing the hot pixel detection 
as soon as it exceeds a given threshold (see Fig. [T]). The latter 
was fixed at four standard deviations of the set of the mean stan- 
dard errors obtained for each position of the eight-point sliding 
window. This method allows us to detect approximately 60-70% 



3 



A. S. Bonomo et al. 2012: Detection of small-size planetary candidates with CoRoT data 



of the most evident hot pixels. After detecting a discontinuity 
in the light curve, we prefer to remove data points in its vicin- 
ity. Indeed, after a simple adjustment of the light curve without 
any clearing of data points in the vicinity of a hot pixel, high- 
frequency variations could still remain and hinder the transit de- 
tection. 

Figure[T]shows one of the light curves selected for our blind 
test that is affected by a clear jump (top panel) and other less 
evident discontinuities. Into this light curve of a r' - 12.8 
G dwarf, namely LRc01_El_0331, transits of a super-Earth with 
Rp = 2.14 /?e and P - 9.8 days were inserted. The red arrow 
indicates the discontinuity identified by our pipeline. In the mid- 
dle panel, the high- and low -pass filtered light curve (3.1 b-c-d) 
is displayed. An artificial dip followed by a jump is seen at the 
time of the hot pixel as a result of the filtering process. Other 
artificial flux drops caused by small undetected discontinuities 
are also discernible. Even though the low-pass filter reduces the 
RMS of the light curve, the individual planetary transits are not 
visible in the light curve. The artificial dip produced by the hot 
pixel identified by our pipeline is corrected after step |3.1[ e (bot- 
tom panel). Despite several artificial dips, in this case the transit 
signal is easily detected: the highest peak at the orbital frequency 
and those corresponding to its harmonics in the BLS spectrum 
are shown in Fig. [2] 

f) Transits are searched for by means of the BLS algorithm 
with the directional correction (cf. Tingley 2003 1. The period 
search was carried out from Pmin - 0.4 up to f max - 30 days. 
The fractionary length of the transit was varied between 0.006 
and 0.09, adopting a number of phase bins «bin,s = 240 ( [Kovacs 
et al.|2002| . The frequency sampling 6v was optimised accord- 



ing to the criterion given by Schwarzenberg-Czerny & Beaulieu^ 
( |2006| l: 6v - l/(Pmax ■ "bins)- This criterion ensures that the sim- 
ulated planets do not remain undetected because of an inappro- 
priate frequency sampling. 




2700 2720 2740 2760 2780 2800 2820 2840 
Time (JD - 2451545.) 



Fig.l. Top panel: the CoRoT light curve LRc01_El_0331, fil- 
tered from cosmic rays (Sect. [3T| a), containing simulated tran- 
sits of a super-Earth with an orbital period of 9.8 days. The red 
arrow indicates the jump produced by a hot pixel and coiTectly 
identified by our pipeline. Middle panel: high- and low-pass fil- 
tered light curve after steps b), c), and d). Planetary transits by 
the super-Earth are not visible to the naked eye. Other artificial 
dips are produced by undetected small discontinuities. Bottom 
panel: light curve after removal of data points in the vicinity of 
the identified discontinuity (Sect. 3.1 e). 
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Fig. 2. Top panel: BLS spectiTim of the filtered light curve 
LRc01_El_0331, containing simulated transits by a super-Earth 
with orbital period P - 9.8 days (see bottom panel of Fig. [TJ. 
The highest peak corresponds to the inserted transit signal, 
which is thus easily detected. Its harmonics are also clearly vis- 
ible. Bottom panel: phase-folded light curve at the period and 
epoch found by the BLS (Sect. [TT 



g) The light curve phase-folded at the period and epoch 
found by the BLS (Fig. [2] bottom panel) is visually inspected 
if events with signal detection efficiency greater or equal than 6 
have been detected. Events below this threshold are generated 
by pure noise with very high probability (Kovacs et al. 2002i. 



In addition, a visual inspection of the original Ught curve fi 
tered from outliers (Sect 



3.1 a) is fundamentally important for 



checking whether a given transit-like feature could be caused by 
undetected hot pixels. The visual inspection checks that there are 
no hot pixels at any epochs of the detected transits or, if any are 
found, that after their removal, the planetary transit is still de- 
tectable and visible in the folded light curve. A complementary 
method is checking possible depth differences when folding the 
filtered light curve at twice the period found by the BLS. Indeed, 
significant depth differences can arise when the transit-like fea- 
ture is produced by a discontinuity in the light curve. 

h) The filtering of an instrumental effect due to the pertur- 
bation associated with the Earth's rotation at the frequency of 
one day"' and its harmonics is carried out. This effect is caused 
by variations of the Earth's infrared emission that produce tem- 
perature changes, even though the way in which these variations 
affect the onboard instrumentation is not well understood yet. 
Such a perturbation gives rise to quasi-sinusoidal signals with a 
period of 12, 24, or 36 h, affecting all CoRoT light curves, al- 
though not in the same way: there is probably a dependence on 
the position of the star on the CCD. 

Mostly, the Earth's perturbation can be neglected and does 
not affect the transit detection. In a few cases, however, it does 
have an effect and may lead to a missed detection. To show this, 
we display in the top panel of Fig. |3] the BLS spectrum of the 
light curve of the r' = 15.0 star LRa01Jil^996, in which we 
simulated transits of a super-Earth with radius - 1.56 R^ 
and period P = 2.1 days. We see that the peak related to the 12 h 
instrumental effect exceeds the power of the transit signal, which 
impedes the transit detection. Indeed, from the visual inspection 
of the phase-folded light curve (Sect. 3.1 g), only a sinusoidal 
variation with a peak-to-peak amplitude of ~ 6 ■ 10 appears 
(see bottom panel of Fig.|3]l. To filter this out, we smoothed the 
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Fig. 3. Top panel: BLS spectrum of the filtered light curve 
LRa01Jil_2996, containing simulated transits by a super-Earth 
with orbital period P = 2.1 days. The 12 h instrumental ef- 
fect overcomes the planetary signal, which impedes its detec- 
tion. Bottom panel: phase-folded light curve at the period and 
epoch found by the BLS, displaying the quasi-sinusoidal vari- 
ation caused by the Earth's perturbation (see text for more de- 
tails). 




tive to discontinuities and/or jumps in the light curves. Indeed, 
more sophisticated filters, such as harmonic fitting (Moutou et alj 
2005; Bono mo & Lanza|2008| ) or wavelet filters (Jenkins 2002]i7 
give rise to deeper artificial flux drops and oscillations in the 
vicinity of an undetected hot pixel. These flux drops, once again, 
can be erroneously identified as transits by the BLS. 



3.3. Optimising fhe extension oftirte median filter window 



Frequency (c/d) 



In the same way as in |Bonomo et al. (p009 ), to find the optimal 
window extension of the sliding median filter, we performed the 
transit search on the light curves of a subset of 500 light curves 
of our blind test by using five different extensions: 0.3, 0.5, 0.75, 
1.0, and 1.25 days. We found that the maximum detection rate 
was reached with a window extension of 0.5 days, which ap- 
parently contradicts the results of [Bonomo et al.| ( |2009| l. These 
authors showed with extensive simulations of light curves with 
solar-like variability and planetary transits that one should use 
the longest possible window extension allowed by the stellar 
activity level to maximise the transit detection. However, their 
light curves were not affected by the instrumental effects that 
are present in CoRoT data. In other words, the short extension 
of 0.5 days allows a better filtering of the CoRoT light curves, 
which are usually affected by sudden jumps and jitter noise 
(Sect. 3.1.e). The longer the median filter window, the stronger 
the filtered light curves are affected by residual instrumental ef- 
fects that hamper the transit detection. 

One may argue that, on the other hand, for long transits, 
a strong filtering with such a short extension could reduce the 
transit signal. However, longer windows did not increase the de- 
tection rate in our exercise as we were searching for small-size 
planets around solar-like stars in short-period orbits, with transit 
durations of typically less than 5-6 h. 
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Fig. 4. Top panel: the same as Fig. [s] after filtering the 12 h in- 
strumental signal (see Sect. |3.1| h). The planetary signal is now 
correctly detected as the most significant peak. Bottom panel: 
phase-folded light curve at the period and epoch found by the 
BLS, distinctly showing the transit signature. 



quasi-sinusoid in phase to determine its amplitude, period, and 
initial phase, and subtracted it from the initial temporal series 
processed up to step e). After which, carrying out again the BLS 
allows us to detect the transit signal because this signal becomes 
the dominant peak in the BLS spectrum (Fig.|4]). 

3.2. Wily use a simple median filter? 

Our choice of using a simple median filter to remove stellar 
variability is motivated by two reasons: the first is the work by 
|Bonomo e t al. (2009 ), who showed that iterative median filters 
( Aigrain & Irwin 2004 ), with an appropriate choice of the win- 
dow extension, perform very well for detections of planetary 
transits; the second is that simple median filters are less sensi- 



3.4. The CoRoT alarm mode 

The CoRoT alarm mode is a CNES operational task that is oper- 
ated at the LAM. It is responsible for detecting transits in CoRoT 
light curves while observations of a given stellar field are ongo- 
ing and uses the so-called Nl level data (Surace et al. 2008| l. 
As soon as transits in a light curve are detected by the alarm 
mode, the temporal sampling of that light curve is changed from 
512 s to 32 s. This is very important for the transit modelling, 
which requires the transit ingress and egress to be well sam- 
pled (e.g., Kipping|2010|), th e study of transit timing variations 
(e.g., Csizmadia et al. |2010 1, and the detection of the secondary 
eclipse (e.g., Alonso et al.|2009 i. Additionally, the alarm mode 
can trigger the ground-based follow-up of the most interesting 
planetary candidates as soon as possible, which is required to 
unveil their nature and, if they turn out to be planets, derive 
the orbital parameters. If a planetary candidate is not detected 
by the alarm mode, it can be discovered by the other detection 
algorithms of the CoRoT team that search for transits on the 
more precise N2 level data (Cab rera et al. |[2009] l. However, in 
this case, the ground-based follow-up has to wait for about six 
months before the target becomes observable again. 

With the exception of CoRoT-21b, all CoRoT planets, in- 
cluding the super-Earth CoRoT-7b, were detected by the alarm 
mode. A replacement of the previous variability filter, a gauging 
filter developed by |Guis & Barge ' (2005 ), with the one described 
in Sect. 3.1 a-d, allowed us to increase the number of planetary 
candidates found by the alarm mode by ~ 20 - 30%. This is 
mainly because the gauging filter is very sensitive to disconti- 
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Table 1. Percentage of detected simulated planets with orbital period 0.7 < P < 20 days for different values of planetary radius. 



Planetary radius 
[Re] 


Detection rate for G and K 
dwarfs with r' < 14.0 [%] 


Detection rate for G and K 
dwarfs with r' < 15.5 [%] 


1.3-2.0 


10.7 ± 1.0 


3.1 ±0.2 


2.0-3.0 


41.7 ± 1.5 


16.8 ±0.6 


3.0-4.0 


76.0+1.6 


44.2 ±0.8 


4.0-5.0 


89.7 ± 1.2 


66.4 ± 0.8 


2.0-4.0 


58.9 ± 1.3 


30.6 ±0.5 
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Fig. 5. Distribution of the simulated planets in the r' magnitude vs orbital period diagram. Green and black dots indicate detected 
and missed planets, respectively. Left: simulations with 1.3 < Rp < 2.5 Rg). Right: simulations with 2.5 < R^ < 5.0 R^. See text for 
explanation. 



nuities caused by hot pixels and shows that light curve filtering 
is crucial for detecting planetary transits, in particular the small- 
size ones. CoRoT-21b was missed by the previous version of the 
alarm mode pipeline that used the gauging filter The replace- 
ment of the variability filter took place in February 2009 during 
the long run LRa02. 

The current alarm mode pipeline is a simplified version 
of the one described in Sect. |3.1[ In particular, it performs 
neither the automatic correction of light curve discontinuities 
(Sect. |3.l]e) nor the filtering of the instrumental effect mentioned 
in Sect. 3.1 h. Indeed, the task of the alarm mode is, in principle. 



the detection of naked-eye transits, hence the complete LAM de- 
tection pipeline (Sect. |3.1[ ) is used only to search for small-size 
planets with N2 data. Nevertheless, even with this simplified de- 
tection pipeline, the alarm mode enables detecting many small- 
size planetary candidates with transits embedded in the noise. 
As an example, it led to the detection of the Neptune CoRoT- 
24c ( |Alonso et al.|[submitted[ l, with a transit S /N lower than 
CoRoT-7b by about a factor of three since it orbits a faint star 
with r' = 15.11 in almost 12 days. 



4. Results and discussion 

4.1. CoRoT detection rate of super-Earths and Neptunes 

With our LAM transit detection pipeline (Sect. [3T[ ), we searched 
for the simulated transits blindly, that is, without knowing in 



advance into which light curves the transits had been inserted 
(Sect. |2|. Figure |5] shows our results in terms of detected plan- 
ets in the r' magnitude vs orbital period diagram. The green and 
black dots indicate the detected and undetected events, respec- 
tively. We considered an event as "detected" when the BLS al- 
gorithm correctly found the period of the simulated transits and 
the signal detection efficiency was > 6. In the left panel of Fig.|5] 
we clearly notice that only few planets with radius Rp < 2.5 R@, 
i.e. 6%, were correctly found. The vast majority of them or- 
bit relatively bright stars and/or have short orbital periods, i.e. 
f < 5 days. In contrast, the right panel shows that about 50% 
of the simulated planets with larger radius, 2.5 < Rp < 5.0 R^, 
were correctly detected. In the latter case, the majority of the 
missed events are the simulated transits with relatively long or- 
bital periods, P > 13 days, which were inserted in faint stars 
with r' > 14.5 (see Fig.|5] right panel). 

Table [T reports the percentage of detections of the simulated 
transits with orbital period 0.7 < P < 20 days for different inter- 
vals of planetary radii. The associated errors were estimated with 
binomial statistics. Below 2 R^, the detection rate ^coRoT is very 
low, ^coRoT ~ 3% by taking all the stars of our blind tests into 
account, and ^coRoT ~ 1 1 % by considering only those brighter 
than r' - 14.0. For the biggest planets with 4.0 < Rp < 5.0 R^,, 
^CoRoT ~ 2/3 (9/10) for r' < 15.5 (< 14.0). By taking small- 
size planets with radii between 2 and 4 P® into account, which is 
the smallest range considered by [Howard et al. (2012 1 to deter- 
mine the planet occurrence rate, our detection rate is 30.6 + 0.5% 
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and 58.9 + 1.3% for stars brighter than r' - 15.5 and 14.0, re- 
spectively. For this range of planetary radius. Table |2] lists the 
fraction of detected planets in eight bins of orbital period, along 
with their own uncertainties estimated, once again, with bino- 
mial statistics. 



Table 2. Percentage of detected simulated planets with radius 
2 < /?p < 4 as a function of orbital period. 



Orbital period 


Detection rate for G and K 


[days] 


dwarfs with r' < 15.5 [%] 


0.7-2.5 


65 ±2 


2.5-5.0 


46 ±2 


5.0-7.5 


36 ±2 


7.5-10.0 


30 ± 1 


10.0-12.5 


24 ± 1 


12.5-15.0 


22 ± 1 


15.0-17.5 


18 ± 1 


17.5-20.0 


15 ± 1 



For each simulated transiting planet we computed a posteri- 
ori the signal-to-noise ratio asS /N = {d/cr) ■ yfn^-, where 6 is the 
transit depth, Wtr the number of data points inside the transit bot- 
tom part, and cr the total noise of the 5 12 s sampled hght curves. 
The latter is estimated as cr = -\/o"w + o"r > where cr^ is the point- 
to-point scatter of the out-of-transit light curve and cTj. its corre- 
lated noise (e.g., Pont et al.|2006 l. cr„ was evaluated in a robust 
way as cr^ - 1.4826 ■ MAD, where MAD stands for median 
absolute deviation, to exclude possible outliers, after removing 
low-frequency variations on time scales longer than 12 hours, cti- 
was estimated by comparing the expected white noise at seven 
different time scales, from 25 min (3 data points) up to 2. 1 h 
(15 data points), with the standard deviation of the sliding aver- 
age over the same time scales (evaluated, once again, in a robust 
way). Since the white noise decreases with the square root of 
the number of measurements, the maximum of the quadrature 
differences between the expected white noise and the measured 
noise, computed for the seven temporal windows, was set to ctj. 
Fig. |6]is an example of the estimated noise in a typical CoRoT 
light curve used in our blind test. In agreement with estimates 
of the correlated noise in CoRoT light curves by |Aigrain et al. 
( 2009| l, we found that its contribution is in general low, ~ 1/4 of 
the total noise. 

Figure |7] displays the radius of the simulated planets in ter- 
restrial units as a function of the signal-to-noise ratio of their 
transits. As in Fig.|5] green dots indicate the correctly identified 
planets while the black dots the missed ones. Figure |7] shows 
that the vast majority of detected transits have > 8 and that 
CoRoT is able to detect transiting small-size planets with an S /N 
much lower than CoRoT-7b, i.e. ~ 1/4 ■ 5/A^coRoT-7b- 



4.2. Obstacles in detecting small-size planets in CoRoTdata 

The vast majority of undetected planets in our blind test were 
missed because their transits have a low S /N. However, some 
events with a relatively high S /N, which should have been quite 
easily identified, were also missed. This motivated us to inves- 
tigate additional obstacles that may contribute to a missed de- 
tection. For this reason, we carefully looked at the few events 
with 5 /A^ > 15 that went undetected (see Fig.|7]i. They represent 
only 5% of the total number of simulated planets with 5 /A^ > 15 
because 95% of them were correctly found. 
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> 0.0003 



0.0002 



0.0001 



0.0000 



5 10 
Number of measurements 



Fig. 6. Measured noise (square) and expected white noise (dia- 
mond) of a typical CoRoT light curve used in our blind test as 
a function of the number of photometric measurements sampled 
at 512 s. In this case, the highest difference is found at 3 mea- 
surements, as indicated by the arrow. 



The first obstacle has been mentioned already (Sect. |3.1| e, 
h) and consists of light curve discontinuities mainly caused by 
hot pixels or pointing displacements that were not identified by 
our pipeline because of their low amplitude. Another obstacle 
is the stellar variability in the case of fast rotators with rotation 
period Piot < 2 - 3 days. In this case, the variability is much 
stronger than the transit signal. Even a very strong filtering is 
unable to reveal the transits because it simultaneously reduces 
their depth, which prevents us from detecting them. We found 
that ~ 70% of the missed transits with S /N > 15 went unde- 
tected by our pipeline owing to unidentified, hence uncorrected 
light curve discontinuities, and ~ 30% because of stellar vari- 
ability. 

While such obstacles are not a major concern for events with 
S /N > 15 since, as already mentioned, only 5% of such plan- 
ets were not detected, they certainly have a stronger impact on 
signals near the detection threshold, i.e. S /N ~ 10. Therefore, 
we suspect they are very likely responsible for the 20% of small 
planets with S /N > 10 that were not found by our pipeline. 

Finally, we recall that the instrumental signal caused by 
Earth's perturbation on 12 or 24 h time scale can in some cases 
also impede the detection of small-size planets, if it is not ad- 
equately filtered. However, this filtering is implemented in our 
automatic pipeline (Sect. [3T| h). 

4.3. Comparison between discovered and expected CoRoT 
planets with /?p = 2 - 4 7?®. 

After deriving the CoRoT detection rate of small-size planets 
from our blind test, we are now able to estimate the expected 
number of CoRoT planets with radius between 2 and 4 R(g from 
the statistics provided by the distribution of the Kepler plane- 
tary candidates. Assuming that these candidates are real plan- 
ets (!Morton & Johnson] [lOTT] !, [Howard et "al~| ([20T2ji found, 
by studying their occurrence rate around G and K dwarfs after 
correcting for the alignment probability, that "the planet occur- 
rence increases substantially with decreasing planet radius and 
increasing orbital period" (see their Figures 5 and 6). They em- 
pirically determined the planet occurrence as a function of plane- 
tary radius and orbital period up to 50 days (cf. their Equations 4 
and 8). In the remainder of the paper we consider only planets 
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Fig. 7. Radius of the simulated planets as a function of the S /N of their transits. As in Fig. |5j green and black dots indicate the 
identified and missed planets, respectively. The vast majority of the discovered transits have S /N > ^ (red dashed line). Note that 
the x-scale in the left and right panels is different. See text for discussion about the transits with S /N > 15 that should be easily 
detectable but were missed in our blind test. 



with 2 < Rp < 4 and 1.2 < P < 17.0 days to choose the 
same bins of planetary radius and orbital period as [Howard et al. 
[ (2012i (see their Figures 4 and 7). Integrating between 1.2 and 
17.0 days gives an occurrence rate of 5.5 + 0.5%. 

The expected number of CoRoT planets with radii of 2-4 7?® 
and orbital periods up to 17 days can be estimated as 



pi 7.0 
2 Jl.2 



dRdP , 



(1) 



where n+[Gv-KV] is the number of G and K main-sequence stars 
with r' < 15.5 in the six CoRoT long runs used in our blind 
test (Sect. [2]i, i.e. those lasting more than 110 days, hence more 
suited to the search for small-size planets; f(R,P) is the fre- 
quency of planets as a function of planetary radius and or- 
bital period determined by |Howard et al. 1 12012); paiP) - 
[(47r2/;3)/(GM*f2)]i/3 fjjg probabihty assuming for 

simplicity = R^, M* = Mq, and circular orbits; ^coRot(^,^') 
is the CoRoT detection rate, which was derived from our blind 
test and is a function of planetary radius and orbital period. The 

r4 nil 

integral Jj ^ f{R^ P) Ptr(P) dR dP represents the fraction of 
transiting Neptune-size planets with P < 17 days according to 
the Kepler planet occurrence (we recall that 20 days is the max- 
imum allowed period in our blind test). The integral in Eq. [T] 
by containing also the factor ^coRot(^, P), gives the percentage 
of such planets detectable by CoRoT. Multiplying the integral by 
n*[Gv-KV] simply allows one to estimate the number of Neptune- 
size planets that CoRoT should have discovered («coRoT[2-4 Re]^ 
see Eq.[T]l. This was evaluated by computing the integrand func- 
tion and multiplying it by n*[Gv-KV] for each of the "cells" that 
range from 2 to 4 /?e and from 1.2 to 17.0 days in Fig. 4 of 
[Howard et al. | ( |2012| l (see Table [3|. The values calculated for 
each cell were then added. 

Because G and K dwarfs are underestimated in the Exodat 
photometric classification, as pointed out by [Gazzano et al.[ 



(2010 1, n^,[Gv-KV] was increased by a factor 1.7^_ajj2V;_KV]_^ 
9516 X 1.7 = 16194 (see Sect. |2]and Fig. 8 in [Gazzano "etaL 
20101. Moreover, a conservative arbitrary error of 30% on this 
quantity was considered because the work by Gazzano et al. 



(2010 1 inevitably relies on a limited number of stars, ~ 1 000. In 



this way, the number of CoRoT planets »coRoTr 2-4 r,,] expected 
from the planet occurrence rate determined by Howard et al. 
([20T2| was estimated to be 12 + 2 (Eq. [T]). The error was evalu- 



ated with a Monte Carlo method and is related to the uncertain- 
ties on «*[Gv-KV], fiR, PX and ^coRot(^, P)- The error on f(R, P) 
was computed with binomial statistics in the same way as de- 
scribed in , Howard et al. [ ( [2012| l, Sect. 3 (see Table|3]l. 

We point out that our estimate of «coRoT[2-4 Re] might be 
slightly underestimated because in our simulations the stellar ra- 
dius and mass were fixed to those of a G2V or K2V star for 
the G and K dwarfs, respectivel}]^ G stars of later spectral types 
would give rise to deeper, hence more easily detectable transits. 
Therefore, it is likely that our detection rate ^coRot(^, P) is mod- 
erately underestimated. However, this effect is partly compen- 
sated for by the slightly lower transit probability for such stars. 

After estimating ncoRoT[2-4 Re]' a comparison with the 
CoRoT Neptunes and Neptune-size planetary candidates is in- 
teresting. For that purpose, the N2 light curves of the six obser- 
vational runs considered for our blind test were analysed with the 
LAM pipeline (Sect. 



3.1 



All small-size candidates previously 
detected by the alarm mode (Sect. [3.4| i and the other six detec- 
tion algorithms used by the CoRoT team, including those dis- 



covered by Ofir et al. (2010 1, were correctly found. The CoRoT 
Neptune-size planetary candidates, regardless of their priority 
and stellar classification, constitute a substantial fraction of the 
total number of candidates (Deleuil et al., in preparation; see, 
e.g.. Fig. 14 in jCarone et al. |2012[ ). 



No information on stellar radii is available in the Exodat catalogue 
and, in the absence of such an information, we preferred to be conser- 
vative. 



8 



A. S. Bonomo et al. 2012: Detection of small-size planetary candidates with CoRoT data 

and number of expected CoRoT Neptunes 



Table 3. Kepler planet frequency f(R,P), CoRoT detection ^coRotC^.^") rate 
«CoRoT[2-4 Re] for each cell of Rp and P as in Fig. 4 of Howard et al. ( 2012 1. 



Planetary radius 


Orbital period 


f(R,P) 




'JCoRoT[2-4 Rj] 


[Re] 


[days] 


[%] 


[%] 




2.0-2.8 


1.2-2.0 


0.037 ± 0.021 


44 ± 6 


0.45 ± 0.31 


2.0-2.8 


2.0-3.4 


0.11 ±0.04 


28 ±4 


0.61 ±0.31 


2.8-4.0 


2.0-3.4 


0.035 ± 0.020 


71 ± 3 


0.49 ± 0.33 


2.0-2.8 


3.4-5.9 


0.51 ±0.11 


19± 1 


1.31 ±0.50 


2.8-4.0 


3.4-5.9 


0.18 ±0.06 


60 ±2 


1.48 ±0.68 


2.0-2.8 


5.9-10.0 


1.0 ± 0.2 


14 ± 1 


1.34 ±0.50 


2.8-4.0 


5.9-10.0 


0.62 ±0.13 


44 ±2 


2.62 ± 0.98 


2.0-2.8 


10.0-17.0 


1.9 ±0.3 


7± 1 


0.85 ± 0.30 


2.8-4.0 


10.0-17.0 


1.1 ± 0.2 


31 ± 1 


2.43 ± 0.83 



Since the statistics provided by [Howard et al. ( 2012)1 does 
not take planets smaller than 2 into account because of 
the current incompleteness of Kepler candidates, we excluded 
CoRoT-7b from this comparison. By considering only the range 
between 2 and 4 R@, we discarded also CoRoT-8b (Borde et 
[al.„2010 ) and CoRoT-24c ( Alonso et al.^submitted ) because all 
of them have radii larger than 4 R(g. The only CoRoT vali- 
dated planet in this range is CoRoT-24b, which has a radius 
of 3.7 + 0.4 R^ and an orbital period of 5.11 days ( [Alonso et 
[aT submitted). We also checked the CoRoT Neptune-size can- 
didates that were not identified as false positives from the light 
curve analysis and follow-up ground-based observations ( Deeg 
[et al.[2009)[Moutou et al.|2009| l. The nature of these candidates is 
unknown, however, because no radial-velocity variation was de- 
tected with the HARPS spectrograph given that they orbit faint 
stars (r' > 13). Therefore, they could be either real planets, in 
which case only an upper Umit can be put on their mass, or 
blended eclipsing binaries not identifiable from the CoRoT light 
curves. For uniformity with our study, we selected the Neptune- 
size candidates orbiting stars with r' < 15.5 that are classified 
as G or K dwarfs by Exodat, in the long runs LRaOl, LRa02, 
LRa03, LRa04, LRcOl, and LRc02. We ended up with three can- 
didates and also added two more candidates whose stars were 
misclassified as giants by the Exodat photometric classification 
but were recognised to be G8V and K3V stars, respectively, 
thanks to low-resolution reconnaissance spectroscopy. However, 
we point out that the transits of two of them are V-shaped, which 
would indicate they are more likely blends. 

In the most favourable case that all these candidates are real 
planets, which cannot be stated with certainty as they have not 
been properly validated, there would be six CoRoT Neptunes. 
This therefore disagrees with the expectation from the Kepler 
mission, i.e. ncoRoT[2-4 Re] = 12 ± 2, at 3 cr. The disagree- 
ment would obviously get worse if a significant fraction of the 
Neptune candidates were actually blends, up to 5 o" for a blend 
fraction of 100%. 

Such a disagreement may be due to several reasons. First 
of all, the uncertainty on f{R, P) might be underestimated be- 
cause the errors on the size of the Kepler planetary candidates 
were not taken into consideration by 'Howard et al. ( 2012| . Such 
errors are related to the uncertainties on the stellar radii in the 
KIC catalogue which are typically ~ 30% (Br own et al. |20lT] l. 
Secondly, the false-positive probability (FPP) of Kepler small- 
size planetary candidates in single systems might be higher than 
predicted by Mort on & Johnson] ( [201 1 ). This has been suggested 
by Colon et al. (2012 1, who foUowed-up four Kepler candidates 



positives ( Colon et al. 



bigger sample of Kep^ 



2012 1. However, Spitzer observations of a 
er candidates would indicate, on the con- 



trary, a l ow FPP ([Desert et al |2012| l. An FPP hi gher than esti 



mated by Morton & Johnsorif l 201 l| l may have led Howard et al. 
( 2012') to overestimate the planet occurrence of small-size plan- 
ets, although at least the Kepler planetary candidates in multiple 



systems are in all likelihood real planets (Latham et al. 2011 



Lissauer et al. 2012| ). An alternative reason might be that a wrong 
value was adopted for n*[Gv-KV] but the latter should be suffi- 
ciently accurate after correcting for the systematic effects of the 
photometric stellar classification and adding an arbitrary error of 
30%. Lastly, the disagreement could also be related to the differ- 
ent stellar populations observed by CoRoT and Kepler In any 
case, this disagreement needs to be investigated in more detail 
in subsequent papers as its explanation goes beyond the scope of 
the present work. 

Despite the above-mentioned disagreement and the fact that 
the exact number of CoRoT Neptunes cannot be precisely de- 
termined for the moment, we point out that their noticeable de- 
ficiency seems to agree with the general trend found in Kepler 
data, i.e. that the frequency of small-size planets increases with 
increasing orbital periods and decreasing planet radii (Howard 



et al. |2012) . If transiting planets with 2.0 < Rp < 4.0 Rfs were 



with radii smaller than 5 R^ and orbital periods less than 6 days 
with multicolour photometry. Half of them turned out to be false 



numerous at short-orbital period, a significant fraction of them 
would have been detected by CoRoT, up to 59% (31%) for G 
and K dwarfs brighter than r' = 14.0 (15.5) (Table[T]i. 



5. Summary and prospects 

We performed a blind test that consisted in simulating transits 
of small-size planets, super-Earths and Neptunes, in real CoRoT 
light curves of G and K dwarfs, and searching for them blindly. 
In this way, we investigated the capability of CoRoT of detect- 
ing small-size planets in short-period orbits. In particular, we 
showed that CoRoT can detect small planets with an S /N as low 
as ~ 1/4 that of CoRoT-7b. We found that the CoRoT detection 
rate is 59% (31%) for planets with radii between 2 and 4 R^ or- 
biting stars with r' < 14.0 (< 15.5) in f < 20 days. The vast 
majority of missed planets went undetected because of a low 
S IN of their transits. However, in some cases, additional instru- 
mental or astrophysical noise may also prevent transits with rel- 
atively high S IN, i.e. S IN > 10, from being revealed. This noise 
is mainly caused by uncorrected light curve discontinuities pro- 
duced by hot pixels and pointing displacements, or short-term 
(< 2 - 3 days) stellar variability. 

By properly taking the CoRoT detection rate for Neptune- 
size planets {2 < Rp < 4 R^) and the transit probability into 
account, CoRoT should have discovered 12 + 2 Neptunes with 
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P < n days around G and K dwarfs of six observational runs, 
according to the Kepler planet occurrence rate (Howard et al. 
2012| . However, between one and six CoRoT Neptunes were 



found in these runs, depending on the fraction of Neptune-size 
candidates that are actually blends. The reason of this disagree- 
ment still needs to be investigated. In any case, by noticing the 
scarcity of close-in Neptune-size planets found by CoRoT, our 
results seem to indirectly support the reality of the trend found 
in Kepler data, i.e. that the frequency of small-size planets in- 
creases with increasing orbital periods and decreasing planet 
radii Poward et al. |2012| l. 

At the end of the CoRoT mission, repeating this com- 
parison by including the next CoRoT long runs and new de- 
tected or validated CoRoT Neptunes will certainly prove instruc- 
tive. Moreover, the radial-velocity follow-up of Kepler small- 
size planetary candidates with the high-resolution spectrograph 
HARPS-N, mounted at the Telescopio Nazionale Galileo (La 
Pahna, Spain), will also unveil their FPP, while determining the 
mass of several of them. The FPP of Kepler candidates is of fun- 
damental importance for studies on the planet occurrence rate 
and could be then compared with the theoretical estimate of 
Morton & Johnson] ( [20n] i adopted by |Howard et al. | ( [20T2] ). 
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